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PLANE ELASTOSTATIC ANALYSIS OF V-NOTCHED PLATES 
by Bernard Gross and Alexander Mendelson 
Lewis Research Center 

SUMMARY 

Solutions are given for several plane elastostatic problems of plates having a v-notch 
on one edge, and subjected to a variety of boundary conditions. The effect of the mag- 
nitude of the v-notch angle and specimen geometry on stress intensity factors Kj and 
Kjj are obtained for unloaded notch surfaces. Examples are given of notch surface dis- 
placements for several selected cases. Analytical solutions are formulated but not 
computed for the case of a loaded notch surface or one with prescribed boundary surface 
displacements. 

Notch opening displacements at the plate edge were measured experimentally and the 
results obtained were in excellent agreement with the computed results. 


INTRODUCTION 

Knowledge of the stress distribution in the neighborhood of a singularity, such as at 
the tip of a v-notch in a loaded plate is of fundamental importance in evaluating the re- 
sistance to quasi-brittle fracture of structural materials. Many different types of spec- 
imens have been used in fracture toughness tests. However, for several reasons con- 
nected with optimization of specimen size and load requirements (ref. 1) single edge 
crack plate specimens loaded in tension or bending are of current interest from the 
standpoint of fracture test method development (ref. 2). 

An analytical solution for a finite width single edge crack specimen loaded in tension 
was given by Gross, Srawley, and Brown (ref. 3), using boundary collocation techniques 
applied to an appropriate stress function derived by Brahtz (ref. 4) and Westergaard 
(ref. 5) and independently by Williams (refs. 6 and 7). Solutions for other finite width 
edge cracked specimens subject to bending or combinations of tension and bending were 
subsequently obtained by Gross and Srawley (refs. 8 to 12). The results are in good 
agreement with those obtained experimentally (refs. 13 to 15 and unpublished data by 
Jones, Bubsey, and Brown of Lewis) and analytically by other methods (refs. 16 to 19), 
as shown in detail in reference 19. 
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(a) Ideal plane crack. (b) Machined notch. 

Figure 1. - In practical specimens the ideal plane crack (a) is simulated by a 
machined notch with a fatigue crack extension (b). The dashed lines 
represent a sharp, finite-angle notch which just includes the machined 
notch and fatigue crack extension. 


The specimens recommended by ASTM Committee E-24 for plane strain fracture 
toughness tests (ref. 2) have machined notches tipped with fatigue cracks (fig. 1(b)) which 
are intended to represent ideal plane cracks (fig. 1(a)). For practical reasons it is de- 
sirable to allow considerable latitude in the form of the machined notch; the only impor- 
tant consideration is that the crack-tip stress field should not be significantly different 
from that of the ideal plane crack which is represented. To establish how much latitude 
in crack -notch configuration can be allowed it is sufficient to investigate the stress fields 
of pointed v-notches as a function of notch angle. Any crack-notch configuration can be 
just bounded by a v-notch of some angle which has its apex at the crack tip (fig. 1(b)). 

The deviation of the stress field of that v-notch from the stress field of a zero-angle 
notch (plane crack) is an upper bound on the deviation of the crack -notch configuration. 
Thus the v-notch angle is useful as a single parameter for characterizing practical 
crack-notch configurations in relation to ideal crack -tip stress fields. The present work 
is concerned with the analysis of v-notched plates as a function of notch angle. 

The first attempt at such an analysis was made in reference 20 for a restricted ge- 
ometry, using the collocation procedures of Gross and Srawley (refs. 3 and 8). The pre- 
sent analysis gives the general solution over a broad range of geometries and in partic- 
ular determines the influence of the notch angle on fracture toughness measurements. 

The elastic stress and displacement solution of a homogeneous isotropic rectangular 
plate of finite dimensions with a v-notch loaded by either symmetric or antisymmetric 
loads is given. The necessary equations for the case of longitudinal shear (mode III) are 
developed but no quantitative solution given. 
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SYMBOLS 


The following list contains the commonly used symbols and their definitions. In 
general, all symbols are defined when introduced. 


A i n 

i,n 


a,b,h, s,w 

E 

G 

K,K r K n 

K c ,K ic’ K nc 

M 


u.u 

r <p 

u , u , u 
x’ y’ z 


A = ri_ + i?„ 
n 'n ’n 

v 


ct . o , a , t 
r’ cp ’ z’ r 0 ’' 


coefficients of the stress function series associated with the n"* 1 eigen- 
value 

specimen dimensions as shown in figs. 2 to 4 
modulus of elasticity 
shear modulus 

stress intensity factor at the notch tip, subscript I, n refers to mode 
of notch extension 

critical value of the stress intensity factor K, Kj, and Kjj 
bending moment 
applied load 

displacements in the radial and tangential directions, respectively 

displacements in x, y, and z directions, respectively 

f v/(l + v) for plane stress conditions 
\ v for plane strain conditions 

■f 

n n complex eigenvalue 
Poissons ratio 
equivalent stress 


T rz ’ V 
V °y’ V T xy 
T yz ’ T zx J 


X 


L components of the stress tensor 


P/bh nominal shear stress at x = 0 for the rectangular double canti- 
lever beam specimen 

stress function 

displacement function 
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ANALYSIS 

The Plane Problem (Modes I and II) 


The solution presented in this paper is for a specimen having a zero radius v-notch 
on one edge. It applies to a homogeneous isotropic material and is in the class of plane 
elastostatic problems of the theory of linear elasticity. 

Plates containing a single edge v-notch are analyzed for symmetric loading (mode I) 
and antisymmetric loading (mode n). The configurations are shown in figures 2 and 3 , 


a 

-LLL I 
B C 



2^- D 
/ a/2 


FRTfl 

o 

Pure tension 
loading 


i 



M 


Pure bending Three point bending 
loading loading (s/w = 4) 



Ct 

h 

Tf 


Double cantilever beam 
specimen. 


Figure 2. - Mode I specimens having a thickness b. 



Figure 3. - Single edge notched plate 
specimen subject to antisymmetric 
loading (mode ID, having a thick- 
ness b. 


respectively. Stress function solutions are obtained in two parts: the first satisfying 
the homogeneous (zero) loading or displacement along the notch surfaces and the second 
satisfying a given load or displacement along the notch surfaces. The most general 
solution is obtained by a combination of these parts. 

We take the origin of coordinates at the tip of the notch as shown in figures 2 and 3. 
Neglecting body forces, the differential equations of the system for both plane stress and 
plane strain reduce to the biharmonic equation 
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( 1 ) 




4 

v X = I 


a 2 , i a , 1 a 2 \/a 2 x, l ax , i a 2 x)_ 
V8r 2 r 3r r 2 i v 2 ) Ur 2 r 3r r 2 a,, 2 ' 


The components of the stress tensor are 


ct - 1 3 X , 1 9X 

r r 2 a 2 r 9r 

r 3 cp 


o = £x 

cp 2 

dr Z 


r + 

r 3r dcp 2 3c p 


It is convenient to define a displacement function \p{ r, cp) such that 


V 2 X= 3^) 

3r 3 cp 


Then the displacements can be written (ref. 7) 


2Gu = - 1* + (1 - K ) r 

r 3r 3 cp 


2Gu = - A 12L + (1 - k) r 2 

^ r dcp 3r 


Assume a stress function x to be given by 

X(r, cp, X) = R e J^Xj(r, cp, X) + x n (r, cp)] 

where R denotes the real part of the bracketed expression, and 


Xj(r, cp,X) = rF(X, cp) 


( 2 ) 


(3) 


(4) 


(5) 
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where 


F(X, cp) ~ £a^ cos Xcp + a-2 cos(x - 2)cp + ag sin Xcp + sin(X - 2)cpJ 


( 6 ) 


and where Xj satisfies identically the conditions of zero notch surface boundary loads 
with X, a^, a 2 , a^, and a^ being complex numbers and Xjj(r, cp) represents a function 
determined from the given nonzero boundary load conditions, if any, along the notch sur- 
faces. 

Assume the displacement function ^ to be 

\jj(r,cp,m) = R e ^j(r,(p,m) + ^ n (r, (p)] (7) 


where 


\lsj(r,(p,m) = r m (d^ cos m cp + d 2 sin m<^) = r m g(<p,m) 


(8) 


and m, dj, and d 2 may be complex numbers. 

Since Xj is independent of Xjj and is independent of ^ it follows from equa- 
tion (3) that 


m = X - 2 


and 


(X - 1) ii- = X 2 F + 

dcp 



(9) 


and 


2 ® '/ / tt ^^11 

V % = r ~ + — 

3r dcp dcp 


From equation (9) we obtain the relation 


d l = 


- 4a 4 
(X - 2) 


and 


d 2 


4a 2 
(X -2) 


( 10 ) 


(ID 
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and it follows that 


A. -2 

^(^*> = ~ £ _ a^ cos(X - 2)<p + &2 sin(X — 2)(^>J + ► 


J 


The homogeneous boundary conditions to be satisfied by the stress function are 


rcp I 




and 


%i r ’f) = %i( r - -f)-° 


From equations (2) and (6) we obtain 


T r Vl = R e 


-r A " 2 (X - 1) II 


dcp 


and 




Satisfying the boundary conditions of equation (13), using equations (14) and (15) 


to 


( 12 ) 


(13) 


(14) 


(15) 

leads 
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and 


-« Ol 

cos X — 
2 

cos (X - 2) — 
2 


A 1 


0 

X sin X — 
2 

(X - 2) sin (X - 2) — 
2 

. 


A 2 


0 


• x a 

sin X — 


a 


sin (X - a) — 
2 


X cos X — (X - 2)cos (X - 2) — 


J 


(16) 


where = 2a^. The resulting eigenequations are 


sin (X n - 1) a = ±(X n - 1) sin a 


(17) 


where 


^n — ^n ^n n ~ 1> 2, 3, . . . , 


are the complex eigenvalues. 

It also follows from equation (16) that 


a 


-cos(X n - 2) _ 


%n 2, n 

cos x - 
n 2 


-sin(X - 2) - 

A 3.n= ^ A 4,„ 

S ‘ nX nI 


The solution for the infinite number of complex roots of these eigenequations was 
obtained by the Newton -Raphson technique and are tabulated in table I for a = 300°. 
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TABLE I. - EIGENVALUES OBTAINED FROM EQUATION (17) FOR a = 300° 




Eigenvalues 

rj + i£ for the even stress function 


n 


5„ 

n 



n 



1 

1. 5122214 

0 

2 

2.4710279 

0. 1418529 

3 

3.6776151 

0.2849014 

4 

4. 8814875 

. 3604963 

5 

6.0840842 

.4136445 

6 

7.2859663 

.4548960 

7 

8.4874029 

.4886817 

8 

9.6885409 

. 5173204 

9 

10. 8894680 

. 5421866 

10 

12.0902400 

. 5641655 

11 

13.2908942 

. 5838618 

12 

14.4914565 

.6017071 

13 

15.6919459 

.6180209 

14 

16. 8923762 

.6330462 

15 

18.0927579 

. 6469723 

16 

19.2930992 

.6599493 

17 

20.4934061 

.6720988 

18 

21.6936841 

.6835201 

19 

22. 8939369 

.6942958 

20 

24.0941684 

.7044951 

21 

25.2943809 

.7141765 

22 

26.4945767 

. 7233902 

23 

27.6947579 

. 7321794 

24 

28. 8949263 

. 7405815 

25 

30.0950830 

. 7486291 

26 

31.2952292 

. 7563511 

27 

32.4953661 

. 7637728 

28 

33.6954947 

. 7709166 

29 

34. 8956156 

. 7778027 

30 

37.2958369 

. 7908716 

31 

37.2958369 

. 7908716 

32 

38.4959383 

.7970851 

33 

39.6960349 : 

. 8031028 

34 

40. 8961263 

. 8089365 

35 

42.0962129 

.8145973 

36 

43.2962956 

. 8200950 

37 

44. 4963746 

. 8254388 

38 

45.6964498 

. 8306372 

39 

46. 8965216 

. 8356977 

40 

48.0965900 

i 

. 8406276 

41 

49.2966561 

.8454334 

42 

50.4967189 

. 8501212 


Tables for other values of ol are given in reference 19. 

The stress function \ can be written separating the even functions and the odd 
functions as follows: 


Xj(r, cp, X) - X IEVEN + Xjqdd 


where 


and 


X IEVEN 


= £ r n [ A l,n cos V + A 2,n cos < X n 
n=l 





> 


X IODD 


w X 

= ^2 r n [ A 3, n sin V + A 4, n sin < X n 
n=l 



Similarly equation (8) is written as 


*I = 


^IEVEN + ^IODD 


where 


and 




IEVEN 


= -4 





A -2 
n 


■a r, cos(A„ - 2) (p 

4 ’ n (* n -2> 


n=l 


( 19 ) 


^IODD " 4 ^ A 
n=l 


V 2 


9 « sin(A -2 )cp 

2."(X -2) ' n 


J 


Since the stresses and displacements are real, they are associated with the real parts of 
the stress function Xj and displacement function \f/j. 

The mode I and n stress and displacement fields are associated with the even and 
odd solutions, respectively. The corresponding stress intensity factors Kj and 
are obtained as follows: 


2 - 7 ], 

K t = \/2n lim r cr (r,0) 
1 r— 0 <P 


L 

= \/2 tj lim r 1 
11 r— 0 


T r <p(r,0) 


( 20 ) 


The stress intensity factors are thus defined to cancel the singularity. As originally used 
by Irwin (ref. 21), the stress intensity factor K is always associated with a zero notch 
angle (crack). This leads to a square root singularity, that is, r = 3/2. When the 
notch angle is no longer zero (crack) but less than 180°, the singularity occurring at the 
notch tip is less than 1/2 (t^ > 3/2, as shown in table VIII). 

Hence 
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Kj = y/2i ?7 1 (?7 1 - 1) 


sin m a sin a 
cos a + 1 I + 1 


1 + cos 77 ^ a 


~\ 


R e A 2, 1 




/sin r].ot sin a \ 

(Vx ~ 2) + Vi[ — cos o') 



1) 

R e A 4, 1 



\ 1 - COS TJjQf / 




- 




( 21 ) 


where rj ^ is the real part of Xj. 

The stress function Xjj is now obtained satisfying V 4 Xjj = 0 and the appropriate 
nonzero notch surface boundary load conditions. For example, if 



( 22 ) 


where kj, kg, k 3 , and k 4 are constants, we can choose the stress function Xjj to be 


X H = r (Ajj + Bjj sin 2 <p + Cjj cos 2 cp + T)^<p) 


And solving for Ajj, Bjj, Cjj, and Djj in order to satisfy equation (22) the stress func- 
tion Xjj is determined, and x = Xj + Xjj satisfies all the field equations and boundary 
conditions . 

It follows from equation (3) and (7) that the displacement function xj / jj must be such 

that 


= 4A n^ + 2D n^ 2 

For the limiting case of a crack ( oi = 360°), the previous equations reduce to those 
obtained in reference 7. 
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B C 

+ 




Figure 4. - Single edge notched plate subjected to longi- 
tudinal shear (mode III). 


THE LONGITUDINAL SHEAR PROBLEM (MODE III) 


For the longitudinal shear problem shown in figure 4, the in plane displacements u r 
and u are assumed to be zero and the longitudinal displacements u z is a function of 
x, y. The notch surfaces are thus constrained to a skew symmetric displacement with 
the surfaces moving parallel to one another. This type of specimen is assumed to be 
very thick so that all planes far enough removed from the free surfaces behave in the 
same way. Superposition of the results of this mode in of loading with that of modes I 
and II is sufficient to describe the most general case of loading on a specimen having a 
v -notch on one edge. 

Introducing a stress function x such that 


rr = - and r = — 

rz cpz 


3 cp 


dr 


the equilibrium equations are identically satisfied and the compatibility equation becomes 


- JLA + 1* + r iLA = v 2 x = 0 


r 3 cp 2 dr 


(23) 


3r 


Equation (23) can be solved by the standard technique of separation of variables, re- 
sulting in 


X = r (Bj cos A cp + B2 sin Xcp) 
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From the boundary conditions along the unloaded notch surfaces, that is, 


T = — = °, cp = ±- 
v z dr 2 


and the fact that x must be an even function of cp, it follows that 


cos — =0, A = (2m - 1) — 
2 m a 


Hence 


X = 


oo 


z 

m=l 


r (2m-l )ir/a 


B 2m-1 


cos(2m - 1) 7 L^£. 

a 


( 24 ) 


where values of m < 1 have been excluded in order that the displacements remain finite. 
While the above solution treats the notch surface as being free (i. e. , x = Xj)> the prob- 
lem of having known shear stresses along the notch surface may be analyzed as well. 

The stress function x is considered to be a combination of two functions X = Xj + Xjj 
where Xj satisfies the homogeneous boundary conditions along the notch surface and Xjj 
satisfies the boundary load conditions along the notch surface. For example, to obtain 
the complete stress function solution assuming a constant shear stress t^ z jj along the 
notch surface, let 


X n = c n r cos <P 

which satisfies the compatibility equation 



= 0 


To determine C-q we have 


3x n 



ar 



cos — = 
2 


T (p zli 
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and X = Xj + Xjj satisfies all the required field equations and boundary conditions. 


BOUNDARY COLLOCATION SOLUTION 

The previous solutions satisfy exactly the boundary conditions along the notch sur- 
faces. To satisfy the boundary conditions along the rest of the boundary, one must find 
appropriate values for the unknown coefficients appearing in the series obtained. This 
can be done by truncating the series as expressed in equation (18) and determining the 
unknown constants such that the stress function and stress function derivative satisfy the 
boundary conditions at a finite number m, of selected boundary stations. In doing this, 
a set of 2m simultaneous algebraic equations is obtained and the first 2m coefficients 
of the stress function series are determined. The truncated stress function series thus 
obtained is considered to be an accurate representation of the actual stress function when 
further increase in the number of boundary stations produce no significant change in the 
first coefficient of the stress function series which represents the dominant term in the 
vicinity of the crack tip. This boundary collocation technique is described in detail in 
references 3, 8 to 12, and 19. 

For included vertex angles a above 300°, the above method worked very well. 
However, as a is decreased below 300° the first coefficient of the stress function os- 
cillated about some value, bracketing this value without actually converging. This dif- 
ficulty was overcome in a manner similar to that used by Hulbert (ref. 22). More bound- 
ary stations were chosen than unknown coefficients leading to an overdetermined system, 
and the resulting set of simultaneous equations was satisfied in the least squares sense. 
All the results reported herein are based on the solution to such overdetermined systems. 
For included angles greater than 300° the results of overdetermined systems were the 
same as for the ordinary system consisting of an equal number of equations and unknowns. 
Those early results referred to in references 3, and 8 to 12 did not use overdetermined 
systems as this was unnecessary. In solving the resulting set of n simultaneous equa- 
tions in n unknowns, a Gauss -Jordon pivotal condensation routine was used. 

Since the most important quantity to be determined is the stress intensity factor K, 
preliminary trials were made to determine the number of equations necessary for eval- 
uating K with sufficient accuracy for the case of pure tension, pure bending, and three- 



TABLE n. - DIMENSIONLESS STRESS INTENSITY 
FACTORS FOR SINGLE -EDGE NOTCHED SPECIMEN 


Included 

vertex, 

angle. 

Notch length to plate width ratio, a/w 

0.2 

0. 3 

0.4 

0. 5 

0.6 

0. 7 

a. 







deg 








r?i-i 

Kjb w 1 







P 



360 

1.085 

1. 614 

2. 369 

3. 539 

5. 537 

9.422 

350 

1.085 

1.613 

2. 369 

3. 541 

5. 538 

9.426 

340 

1.087 

1.618 

2. 374 

3. 549 

5. 549 

9.445 

330 

1.097 

1. 630 

2. 389 

3. 569 

5. 579 

9.491 

300 

1. 169 

1. 724 

2. 520 

3. 756 

5. 859 

9. 979 

270 

1. 366 

1. 987 

2. 888 

4. 297 

6. 736 

11. 515 

240 

1.804 

2. 595 

3. 766 

5.630 

8. 934 

15. 551 


K i 

Pure bendinp-! 

7 

bw 1 







6M 



360 

0. 837 

1. 093 

1.414 

1. 877 

2.629 

4.041 

350 

. 836 

1.093 

1.414 

1. 876 

2. 627 

4.044 

340 

. 839 

1.095 

1. 416 

1. 878 

2.632 

4.048 

330 

. 844 

1. 100 

1. 422 

1. 885 

2.640 

4.062 

300 

. 895 

1. 155 

1.484 

1. 965 

2. 752 

4.250 

270 

1.034 

1. 314 

1.678 

2.218 

3. 139 

4. 873 

240 

1. 346 

1. 680 

2. 146 

2. 860 

4. 112 

6. 532 


j 

Kjb w 

V 1 




i'V'****. 


6P 


360 

0. 777 

1.013 

1. 318 

1. 767 

2. 504 

3. 893 

350 

. 776 

1.012 

1. 319 

1. 766 

2. 504 

3. 898 

340 

. 779 

1.014 

1. 320 

1.767 

2. 505 

3.903 

330 

. 785 

1.021 

1. 325 

1. 775 

2. 517 

3. 922 

300 

. 832 

1.072 

1. 390 

1. 857 

2. 636 

4. 112 

270 

.963 

1.226 

1. 578 

2. 118 

3.012 

4. 736 

240 

1.261 

1. 573 

2.043 

2. 770 

3.984 

6. 395 




TABLE IB. - DIMENSIONLESS STRESS INTENSITY FACTORS FOR AN 
EDGE-NOTCHED RECTANGULAR DOUBLE -CANTILEVER BEAM 
SUBJECT TO SPLITTING FORCES 


Included Maximum beam 
vertex depth to notch 
angle, depth ratio, 
a, h/a 


Notch depth to plate width ratio, a/w 


0.667 0.500 0.400 0.333 0.250 0.222 0.200 


8.320 6.565 6.560 6.485 6.587 

14.829 7.202 5.810 5.798 5.798 

22.209 9.896 6.033 5.695 5.678 

29.623 13.039 6.958 5.913 5.730 

8.344 6.601 6.571 6.593 6.606 

14.837 7.211 5.822 5.815 5.810 

22.213 9.899 6.033 5.699 5.709 

29.630 13.038 6.952 5.910 5.752 

8.360 6.619 6.634 6.611 6.732 

14.853 7.215 5.830 5.812 5.810 

22.240 9.903 6.035 5.696 5.678 

29.589 13.048 6.950 5.906 5.789 

8.383 6.639 6.636 6.641 

14.904 7.230 5.835 5.822 5.822 

22.320 9.928 6.040 5.705 5.689 

29.775 13.085 6.956 5.904 5.787 

15.509 7.426 5.964 5.947 5.961 

23.263 10.253 6.168 5.781 5.767 

31.040 13.541 7.122 6.071 5.881 

26.333 11.300 7.741 6.509 6.707 

35.139 14.985 9.732 7.654 6.307 


point loading. Over the range of geometries analyzed it was indicated that an overde- 
termined system of 52 equations and 40 unknowns were sufficient. For the double canti- 
lever beam antisymmetrically loaded, it was found that a set of 64 equations and 40 un- 
knowns were sufficient. 


Numerical Results 

Calculations were performed for mode I notch opening for the cases of pure tension, 
pure bending, three -point loading, and the double cantilever beam subjected to splitting 
forces as shown in figure 2. Calculations were also performed for the antisymmetric 
mode n notch opening as shown in figure 3. The required boundary conditions for all 
these cases are given in the appendix. 

Preliminary results indicated that for pure tension, pure bending, and three -point 
loading, increasing the height to width ratio h/w beyond 1. 2 (fig. 1) for a > 300° and 
beyond 1. 4 for 300° > a > 240° produced no change in the stress intensity factor. These 
h/w values were therefore used in all the subsequent calculations. 

The values of the computed stress intensity factors Kj and Kjj are given in 
tables II to IV for various values of notch angle a and notch depth ratios a/w. For the 


TABLE IV. - DIMENSIONLESS STRESS INTENSITY FACTORS 
FOR SINGLE-EDGE NOTCHED SPECIMEN SUBJECTED TO 


ANTISYMMETRIC LOADING, MODE II 


Included 

Maximum beam 

Notch depth to plate width ratio, a/w 

vertex 

angle, 

depth to notch 
depth ratio, 

0. 400 

0. 333 

0.286 

0.200 

a, 

deg 

h/a 


K n 

bh 




2 

Pa 



360 

0. 5 

0. 353 

0. 353 

0. 353 

0. 352 


1.0 

.500 

.500 

. 500 

.500 


1. 5 

.614 

.612 

,.612 

.612 

350 

0.5 

0. 186 

0.186 

0. 186 



1.0 

.401 

.401 

.401 



1. 5 

. 543 

. 540 

. 540 


340 

1.0 

0.278 

0.278 

0.278 



1. 5 

.456 

.454 

. 454 


330 

1.0 

0. 124 

0.124 

0.124 



1.5 

. 348 

. 346 

. 346 
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double cantilever beam, the beam depth to notch depth ratio h/a was also varied as 
shown in table HI. 

Typically computed equivalent stress contours and notch displacement values for the 
case of three-point loading are shown in figures 5 and 6. Tables and figures are given 
in reference 20, covering a much wider range of geometries and loading. 


Stress 

ratio, 



stress conditions. 



Location, x/a 

(b) Dimensionless equivalent stress contours for plane 
strain conditions. 

Figure 5. - Edge notched beam subjected to three point 
loading. Included vertex angle, 300°; Poissons ratio, 

0. 30; notch length to plate width ratio, 0. 30; a = 6P/bw. 



Figure 6. - Three-point loading elastic notch edge dis- 
placements for plane stress conditions. Included 
vertex angle, 300°. 


EXPERIMENTAL MEASUREMENT OF DISPLACEMENT 

For comparison with the analytical results, crack surface displacement measure- 
ments were made on several v-notched plate specimens of 7075T651 bare aluminum, 
which were subjected to three -point bending. Measurement points were symmetrically 
located on opposite sides of the edges of the v-notch. Figures 7 to 9 shows the equip- 
ment used. Numbers appearing in figures 8 and 9 refer to the various components in- 
dicated in figure 7. The bend specimens were 9 -inches long by 2 -inches wide by 1-inch 
thick with a span length -to -width ratio of 4 to 1. Nine specimens were tested corres- 
ponding to a values of 330°, 300°, and 270°, and a/w ratios of 0.4, 0.5, and 0.6. 
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Loading unit 
crossarro vertical 
~ displacement 
control and 


| load recorder 


10 


Power supply 1 
supplies excitation 
voltage and bridge 
balan ce 7 


Excitation voltage 
and load cell output 


Excitation voltage 
and displacement 
gage output 


| Crossarm | n 


Standard 
cell output} 


I 


Digital readout f 
microvolt mete r r 

[[Digital printouf 

Standard cell output 
Load cell bridge excita- 
tat ion voltage 
Displacement gage bridge 
excitation voltage 
load cell output 
Displacement gage output 



Displacement 

gage 

1 


5 Calibrated sensing 
compressive load unit 


Power supply 2 
supplies excitation 
[voltage and bridge 
balance 6 


Figure 7. - Schematic of experimental test rig. 



Figure K. - Close up of loading fixture am instrumented specimen. 
(Numbers correspond to those shown m fig, 7.) 


co 









Figure 9. - Loading frame, fixtures, and instrumentation. (Numbers 
correspond to those in fig. 7.) 


The modulus of elasticity of the aluminum was taken as 10.4x10° psi. The equipment 
used and the experimental procedures are described in detail in reference 19. 


RESULTS AND DISCUSSION 

Results are presented in tables 13 to V. The value of the first term of the truncated 
even stress function R e (A2 j) or R e (A^ is directly proportional to the stress inten- 
sity factor Kj or Kjj. Tables II to IV contain the nondimensional stress intensity 
factors for mode I and mode n. 

For all mode I cases considered the value of K increased a maximum of 1 percent 
in going from a = 360° (crack) to 330°. The square root singularity associated with the 
360° crack changed from 0. 5000 to 0.4996 in going from 360° (crack) to 330°. One can 
conclude on examining the results of table HI that as long as a is 330° or greater the 
difference in the stress intensity factor Kj from that for a crack is very small. Table IV 
contains the nondimensional mode II stress intensity factors. Clearly, small changes in 
a produce large changes in K^. 
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TABLE V. - COMPARISON OF EXPERIMENTALLY 


MEASURED AND COMPUTED PLANE STRESS 
DIMENSIONLESS DISPLACEMENT VALUES 
FOR THREE -POINT BENDING 


[Span to width ratio, 4 to 1; cr = P/bw. ] 


Specimen 

number 

Included 

vertex 

angle, 

a, 

deg 

Notch depth 
to plate 
width ratio, 
a/w 

Dimensionless 
plane -stress 
displacements, 
2EUy/a<j 

Experimental 

Computed 

1 

330 

0.4 

8.69 

8.72 

2 



. 5 

11.86 

11. 89 

3 



.6 

17. 71 

17. 86 

4 

3( 

30 

0.4 

8. 90 

8.92 

5 



. 5 

12. 17 

12. 19 

6 



.6 

18.07 

18.37 

7 

2 

70 

0.4 

9.46 

9. 55 

8 



. 5 

13.01 

13. 12 

9 



.6 

19. 80 

20.00 


Figure 5 contains typical closed contour curves of dimensionless equivalent stress. 
One obtains from these contours a semiqualitative look into possible regions of plastic 
flow and their shapes. Straight lines have been drawn between points of constant equiva- 
lent stress along rays emanating from the crack tip at 10° intervals. 

A typical dimensionless plane stress y displacement curve along the notch surface 
is plotted in figure 6 for mode I. Displacements were computed along the notch at inter- 
vals of one -tenth the relative notch depth r/a. Corresponding plane strain displacements 

o 

may be obtained by multiplying the plane stress y displacement by (1 - v ). 

Confidence in the correctness of the results obtained by the method of the present 
analysis was derived from many comparisons, both with the experimental results ob- 
tained herein as well as experimental results of other authors for the crack (o’ = 360°) 
problem. Additional verification was obtained when comparisons were made with solu- 
tions to several special cases by other investigators using different analytical techniques. 

Table VI shows a comparison of the dimensionless stress intensity factor obtained by 
the present method for the cases of pure tension, pure bending, and three-point loading 
with experimental results (refs. 13 to 15) and analytical results (ref. 18) of other inves- 
tigators. Good agreement is obtained. 

Table VTE shows a comparison of the plane stress dimensionless displacement ob- 
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TABLE VI. - COMPARISON OF DIMENSIONLESS STRESS INTENSITY FACTORS FOR SPECIMEN 


SUBJECTED TO PURE TENSION, PURE BENDING, AND TH RE E - POINT LOADING 


[Vertex angle, 360°; Mode I (see fig. 2). ] 


Type of data 

Source 



Notch depth to plate width ratio, a/w 






0. 10 

0. 15 

0.20 

0.25 

0. 30 

0. 35 

0.40 

0.45 

0. 50 




p p p 

Pure tension; dimensionless stress intensity factor, K^wb /P 


Experimental 

Ref. 

15 

0. 556 

0. 816 

1. 180 

1. 735 

2. 571 

3.775 

5.436 

7.641 

10.477 

Experimental 

Ref. 

13 

.65 

1.00 

1.40 

1.97 

2.80 

4.20 

6. 18 

8. 90 

12. 50 

Collocation 

Ref. 

5 

.445 

. 758 

1. 180 

1. 768 

2.603 

3. 813 

5. 596 

8.276 

12. 399 

Analytical 

Ref. 

18 

.443 

. 747 

1. 164 

1. 751 

2. 592 

3. 813 

5.606 

8.284 

12. 363 




2 2 

Pure bending; dimensionless stress intensity factor, K^b w' 

3 /M 2 ; 


Collocation 

Ref. 

8 

12.4 

18. 5 

25. 3 

33.2 

42. 8 

55.2 

71.4 

92. 7 

123.0 

Experimental 

Ref. 

14 

11.8 

17. 4 

24.2 

32. 15 

41.9 

53.9 

68. 6 

88.9 

118.0 

Analytical 

Ref. 

18 

12. 37 

18.28 

24.94 

32.95 

42. 85 

55. 57 

72. 14 

94. 5 

125. 5 



2 

Three-point loading; dimensionless stress intensity factor, Kjb 

2 w 3 /M 2 ; 





span to width ratio, 4/l 







Collocation 

This work 



21. 78 


36. 89 


62. 50 


112. 34 

Experimental 

(a) 




20.26 


37.62 


61. 75 


110. 77 


Unpublished data obtained by Jones, Bubsey, and Brown of Lewis. 


TABLE VIL - COMPARISON OF PLANE-STRESS DIMENSIONLESS DISPLACEMENT ACROSS NOTCH SURFACES 
FOR THE THREE-POINT LOADING SPECIMEN AND RECTANGULAR DOUBLE-CANTILEVER 
BEAM SUBJECT TO SPLITTING FORCES 


[Vertex angle. 360°. see fig. 2. | 


Type of data 

Source 

Gage location 

Plate width 


Notch depth to plate width ratio, a/w 




x/w 

y/w 

to beam 
height 

0.20 

0. 30 

0. 40 

0. 50 

0.60 

0. 70 

0. 80 





ratio. 












w/h 










Three -point loading; dimensionless displacement, 2EbUyW/6M 




Collocation 

This work 

-a/w 

±0 

5/6 

1. 159 

2.075 

3.497 

5. 942 

10. 672 

21. 446 

51. 370 

Collocation 

This work 

-a/w 

±.10 

5/6 

1.166 

2.077 

3.497 

5.943 

10.672 

21.445 

51. 399 

Experimental 

(a) 

-a/w 

±.10 

— 

1. 176 

2. 124 

3. 556 

6.032 

10.29 

21.09 

51.53 

Rectangular double 

-cantilever 

beam subject to splitting forces; dimensionless displacement. 

2Ebu /P 

Collocation 

This work 

-(a/w+l/4) 

±0 

5/3 


23.47 

34. 73 

50. 13 

87. 18 

163. 39 

395. 7 


Experimental 

(a) 

-(a/w+1/4) 

±0 

5/3 

24.21 

35.44 

53. 93 

88. 83 

161. 14 

380.0 



a UnpubIished data obtained by Jones, Bubsey, and Brown of Lewis. 
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TABLE Vffl. - TABLE OF FIRST 


EIGENVALUES FOR MODES I 
AND n FOR VARIOUS 
INCLUDED VERTEX 


ANGLES 


Included vertex 
angle, 
deg 

First eigenvalues, 77 ^ 

Mode I 

Mode II 

360 

1. 500000 

1. 500000 

350 

1. 500053 

1. 529355 

340 

1. 500426 

1. 562007 

330 

1. 501453 

1. 598192 

300 

1.512221 


270 

1. 544484 


240 

1. 615731 



tained by the technique presented herein with that obtained experimentally (unpublished 
data obtained by Jones, Bubsey, and Brown of Lewis) for the cases of three-point bend- 
ing and rectangular double -cantilever beam subject to splitting forces. For the rec- 
tangular double -cantilever case the displacements were extrapolated linearly to the 
specimen edge [x = -(a + w/4), y = 0]. Once again good agreement is obtained. A more 
detailed comparison of the analytical and experimental results is given in reference 19. 

For included vertex angles a less than 360°, no experimental or analytical data 
have heretofore been available. Experiments, as previously described were therefore 
made for a three -point bend specimen, and the results were compared with the analysis 
presented herein for the case of plane stress and are shown in table V. Results are 
given for a values of 330°, 300°, and 270° and notch depth to plate width ratios a/w 
of 0. 4, 0. 5, and 0.6. Excellent agreement is obtained over the range of a/w ratios. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 1, 1970, 

731-25. 
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APPENDIX - STRESS FUNCTION AND STRESS FUNCTION DERIVATIVE 
ALONG THE BOUNDARY OF v-NOTCH SPECIMEN 


For pure tension as shown in figure 2 


Along AB 


Along BC 


Along CD 


x = o = o 

0 X 


x = J_ 

P bw 



= 0 

ay 


X _ w_ b 9x = 1 

P 2b P dx 


For pure bending as shown in figure 2 
Along AB 


X = 0 12U 0 

Sx 


Along BC 


+ ^ + + « + i*=o 

M w 3 y 6 2 2 6/ w 2 y 2 2j 9y 

Along CD 


bx = i lx = o 

M 9x 
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n 


For three -point loading as shown in figure 2 
Along AB 

X= 0 0 

dX 


Along BC 


X_ _ 6 

P bw' 


bw 3 \ 2 / 


3 2 _2 

- — + — (w - 2a) + — (w - a) + 

6 4 2 


a 2 /w a\ 
2 \2 3/ 


b 0x = 

p 0 y w 3 


3 2 

x . x 


- — + — (w - 2a) + — (w - a) + 

64 2 2 V 2 3 ; 


a /w a' 


Along CC’ 


Along C’D 


X = s - 2y 3x _ o 

P 4b 9x 


x = JL 

P 4b 


s - 2y 


(y - y 0 )' 


= 0 
9x 


For the rectangular specimen subject to splitting forces as shown in figure 2 
Along AB 


r 2 

X_ _ P (y + g) sin y_ sin 0 + sin y sin(y + 9 ) 
P b L 2y + sin 2y 

2 

(0 + y) cos y cos 0 - cos y sin(y + 0) 
2y - sin 2y 

+ _e_ / sin 2 9 - 20 cos 2y 
2b \2y cos 2y - sin 2y 
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where 


277 - O' 

4 


e = 


u . a 

h + a tan — 
2 

tan 2 y 


p = V (h - r sin cp)^ 1 + (e + a + r cos cp 

6=y- tan' 1 / h-rsiny \ 

\e + a + r cos cp) 

and the stress function derivatives is obtained from 


* .iX= b Vx . j 
p ax p 

where f is a unit vector normal to boundary. 

Along BC 

x = / x + a\ ax = o 

P \ b / 3y 

Along CD 

x, _ w _b ax = j 

p b p ax 


For the Antisymmetric loading as shown in figure 3 
Along AB 
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Along BC 


, i a 

h + a tan — 

X_ = 2 _b dx = 1 

P 2b P 9y 

Along CD 



For the longitudinal shear (mode HI) as shown in figure 4 the stress function require- 
ments along the specimen boundary is as follows 

Along AB 


X = T xz i sin (p + a tan 


Along BC 


* =T xzl( h + atan fj 


Along CD 


* =T xz2 rsin <P + h ( T xzl ' 
The detailed derivation of these conditions 


t xz2> + T xzl a tan | 
is given in reference 19. 
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